Collinear versus non-collinear magnetic order in Pd atomic clusters 
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We present a thorough theoretical assessment of the stability of non-collinear spin arrangements 
in small palladium clusters. We generally find that ferromagnetic order is always preferred, but that 
antiferromagnetic and non-collinear configurations of different sorts exist and compete for the first 
excited isomers. We also show that the relative stability of all these states is rather insensitive to 
the choice of atomic configuration for the pseudopotential used and to the approximation taken for 
the exchange and correlation potential. This result stands in stark contrast with the situation found 
for the bulk phases of Palladium. 
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The magnetic properties of free-standing atomic clus- 
ters of 3d TM elements have been intensively scrutinized 
during the last two decades. Two different but related 
phenomena have specifically been discussed and essen- 
tially unravelled. The first is the modification of local 
magnetic moments as compared with the values found 
in bulk materials. The second is the competition be- 
tween the possible ferromagnetic, antiferromagnetic and 
non-collinear arrangements of the local spins, as well as 
its interplay with the geometry of the nanostructure. fn 
the case of ferromagnetic elements like Fe, Co and Ni, 
the increase of the average cluster magnetic moment can 
be easily explained in terms of the reduced atomic co- 
ordination in the low-dimensional regime, with oscilla- 
tions associated to structural (symmetry) changes, |lj 
The case of antiferromagnets like Cr and Mn is much 
more complex. Atoms of these elements may display 
large magnetic moments, since they have a large num- 
ber of d-holes susceptible to be polarized. On the other 
hand, clusters of these atoms may display tiny average 
magnetizations due to the tendency of their atomic mo- 
ments to align in antiparallel directions. The structure 
plays also a fundamental role in the magnetic behavior 
of these clusters, since it may originate magnetic frus- 
tration. A conventional example of magnetic frustration 
in a classical spin system appears when atoms positions 
form triangular motifs. The studies of these classical sys- 
tems show that magnetic frustration frequently leads to 
non-collinear configurations of the local spin moments. 
The latest theoretical studies reported in the literature 
show that non-collinear arrangements of quantum spins 
also appear as the ground or as some of the first isomers 
of clusters of 3d atoms, including not only Cr and Mn, 
but also Fe, Co and Ni.0,|3,|j,|3,|g 

All materials made of 4c? TM elements are paramag- 
nets, in contrast to those of the 3d row. A natural ques- 
tion thus arises of whether small clusters of Ad elements 
may show low-lying magnetic states of collinear or even 
non-colinear nature. Bulk palladium, being a paramag- 
net in the brink of becoming a ferromagnet, presents one 



of the most intriguing and controversial magnetic behav- 
iors in nature. [3 It is therefore not surprising that the 
very few experimental and theoretical studies published 
so far try to clarify whether Pd clusters of given sizes are 
magnetic or not, and what is the order of magnitude of 
their average magnetic moment. From the experimental 
side, most of the reports agree that only very small clus- 
ters have a net magnetic moment H H Eolllll , with the 
exception of Shinohara and coworkers, [12| who found 
noticeable magnetic moments at the surface of Pd par- 
ticles as big as 79 A. From the theoretical side, there 
is also consensus that very small Pd clusters are indeed 
magnetic. [H Q III III 111 El Futschek et al. have 
studied recently small Pd clusters using Density Func- 
tional Theory (DFT) in the collinear framework, within 
a fixed-moment mode. They have found that multiple 
spin isomers exist for each cluster size with very small 
energy differences. Interestingly, some of these compet- 
ing isomers present ferromagnetic order, while others dis- 
play antiferromagnetic alignments, with possible frustra- 
tion. Although Pd has tendency to ferromagentic order, 
this fact strongly points out to the possible existence of 
non-collinear magnetic structures, as a mechanism to re- 
lease the frustration and competition between the differ- 
ent magnetic solutions. 

We report in this article a thorough Ab initio study of 
the magnetic behavior of small palladium clusters Pd^, 
with N ranging from 3 to 7. We have performed a si- 
multaneous optimization of the geometric and magnetic 
degrees of freedom fully allowing for non-collinear spin ar- 
rangements. This consists, to the best of our knowledge, 
the first study of non-collinear magnetism in Ad atomic 
clusters. Moreover, a debate currently exists on the ac- 
curacy of the Local Density Approximation (LDA) (2(J 
versus the Generalized Gradient Approximation (GGA) 
[2l| for the determination of the ma gne tic behavior of 
low-dimensional Pd systems [II El El El E3 The 
present letter also assesses the reliability of both approx- 
imations for the case of free-standing Pd atomic clusters. 

We have performed our calculations using the code 
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FIG. 1: Illustration of the ground state structures of the dif- 
ferent clusters here studied and average interatomic distances 
(in A) within GGA1. 

SIESTA. SIESTA is a DFT method that employs 
linear combination of pseudoatomic orbitals as basis 
set. The electronic core is replaced by a nonlocal 
norm-conserving Troullier-Martins |2 5| pseudopotential 
that may include nonlinear core correction terms. The 
code allows to perform, together with the electronic cal- 
culation, structural optimization using a variety of al- 
gorithms. It also allows to simulate non-collinear spin 
arrangements both in the LDA and in the GGA approx- 
imations. |26j 

In the present calculation, we have also used a variety 
of pseudopotentials to test their effect on free-standing 
clusters and their corresponding transferability. We have 
generated three different pseudopotentials using LDA. 
The first (LDAl) was built with the electronic config- 
urations 5s 1 , 5p° and 4c? 9 , and core-corrections matching 
radius r c — 2.00 a.u.; the second (LDA2) was identical to 
LDAl, but with r c = 1.2 a.u.; the third had a closed-shell 
atomic configuration (5s , 5p° and Ad 10 ) and r c — 1.2 a.u. 
We have also generated two GGA pseudopotentials with 
electronic configuration 5s 1 , 5p° and 4d 9 , and r c = 2.0 
or 1.2 a.u. (GGA1 and GGA2, respectively). In all five 
cases, the cutoff radii of the s, p and d orbitals were 
taken at 2.30, 2.46 and 1.67 a.u., respectively. We have 
described valence states by a double-^ polarized basis set 
(e.g.: two different radial functions for s and d orbitals 
and a single one for p orbitals). We have taken an energy 
cutoff of 150 Ry to define the real space grid for numer- 
ical integrations, but we checked that higher cutoffs did 
not alter the results. We have carried out the structural 
optimization using a conjugate gradient algorithm, where 
we have set the tolerance for the forces at 0.003 eV/A, 
with eventual double-checks using 0.001 eV/A. 

We have found that the five pseudopotentials provide 
similar results when applied to an isolated palladium 
atom, being the eigenvalues of the ground state and 
different excited states slightly better reproduced with 
LDAl and GGA1 (both had r c = 2.00 a.u.). However, 
we have observed that they give rise to different magnetic 
behaviors when applied to the bulk fee material. All LDA 
approximations give a lattice constant equal to 3.90 A, 
while all GGA predict it to be equal to 4.01 A. LDAl 
gives a ferromagnetic ground state with M » 0.54/is, 
while LDA2 and LDA3 predict the ground state to be 
paramagnetic. Finally, both GGA pseudopotentials lead 
to a ferromagnetic ground state with M ~ 0.48/zg. These 



TABLE I: Bindig energy of the ferromagnetic clusters in 
meV/atom. 



N 


LDAl 


LDA3 


GGA1 


Ref.flSl 


RefJl9] 


3 


1.755 


1.326 


1.289 


1.203 


1.250 


4 


2.293 


1.942 


1.769 


1.628 


1.675 


5 


2.502 


2.168 


1.933 


1.766 


1.805 


6 


2.721 


2.401 


2.110 


1.919 


1.949 


7 


2.791 


2.452 


2.155 


1 .953 


1 .985 



results highlight the importance of testing all the differ- 
ent pseudopotentials for atomic clusters considered here. 

Notice that we have not kept fixed the magnetic mo- 
ment in our simulations of the Pd^ clusters, but rather 
have allowed it to vary freely during the non-collinear 
iterative selfconsistency process, in contrast to previous 
authors. Moreover, while we can not rule out that we 
may have missed low lying solutions, we have endeav- 
ored to minimize this risk by feeding a large variety of 
non-collinear seeds for each cluster. This effort has al- 
lowed us to find a rich and complex family of metastable 
solutions, that was absent in previous works. We finally 
note that we have repeated all calculations with the pseu- 
dopotentials LDAl, LDA3 and GGA1. 

We have found that all clusters, except Pd6, share the 
same collinear magnetic ground state, with a total spin 
of 2 /j,b, in agreement with previous authors 003. We 
should stress that all the tested pseudopotentials provide 
the same ground state, in stark contrast to the situation 
that arose for the bulk material. Moreover, we have found 
very similar inter-atomic distances for all Pdjy clusters, 
using whichever pseudopotential. These distances also 
agree with those obtained by Kumar and Futschek within 
a range of 1 per cent. The geometry of the ground state 
and the average interatomic distance of the Pdjv clusters 
is displayed in Fig. 1, where we show that these range 
from 2.55 A in Pd3 to 2.71 A for Pdy. We have written 
the binding energies of the different clusters in Table I. 
The table shows that GGA1 gives slightly smaller values 
than LDAl and LDA3, as otherwise expected. Moreover, 
the binding energies predicted by GGAl are very similar 
to those obtained by Kumar, who also used the GGA 
(within an ultrasoft pseudopotentials, plane waves code) 
and by Futschek et al., who used the all-electron VASP 
code, but did not statp thp annmYimatinn employed. 




FIG. 2: Illustration of the non-colinear magnetic solutions 
for Pd 4 NC2 (LDAl), NC2 (LDA3) and NCI (GGAl). The 
arrows are proportional to the size of the atomic moments. 
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TABLE II: Different solutions obtained for the Pdjv clusters. We provide the absolute values of the atomic magnetic moments, 
the total magnetic moment in the cluster p, (both in units of fis) and the excitation energy per atom (in meV). For N=5 and 7, 
the first two values of the atomic moments correspond to the axial sites, whereas the last ones correspond to the planar sites. 





LDA1 






LDA3 






GGA1 








Local moments 




AE 


Local moments 




AE 


Local moments 




AE 


N=3 


Ferro. 


(0.67x3) 


2 





(0.67x3) 


2 





(0.67x3) 


2 





AF 














(0, 0.30, -0.30) 





28 


Radial 














(0.18x3) 





28 


Para. 


(0x3) 





68 


(0x3) 





34 


(0x3) 





75 


N=4 


Ferro. 


(0.50x4) 


2 





(0.50x4) 


2 





(0.50x4) 


2 





NCI 














(0.29,0.29,0.29,0.29) 





12 


NC2 


(0.35,0.24,0.24,0.35) 


0.25 


9 


(0.25,0.28,0.28,0.25) 


0.03 


10 








AF1 


(0.32, 0.32, -0.32, -0.32) 





26 


(0,23, 0.23, -0.23, -0.23) 





30 


(0.29, 0.29, -0.29, -0.29) 





25 


AF2 


(0.41, 0, -0.41, 0) 





40 


(0,31, 0, -0.31, 0) 





31 


(0.38, 0, -0.38, 0) 





36 


Para. 


(0x4) 





86 


(0x4) 





59 


(0x4) 





78 


N=5 


Ferro. 


(0.43,0.43,0.38x3) 


2 





(0.40x5) 


2 





(0.42,0.42,0.39x3) 


2 





AF1 


(0, 0, 0.43, - 0.43, 0) 





22 


(0, 0, 0.33, - 0.33, 0) 





19 


(0, 0, 0.39, - 0.39, 0) 





18 


AF2 


(0, 0, 0.48, -0.24, -0.24) 





27 








(0, 0, 0.44, -0.22, -0.22) 





19 


Radial 


(0, 0, 0.29x3) 





35 








(0, 0, 0.27x3) 





28 


Para. 


(0 x5) 





63 


(0 x5) 





41 


(0 x5) 





55 


N=6 


Ferro. 


(0.33x6) 


2 





(0.33x6) 


2 





(0.33x6) 


2 





Para. 


(0x6) 





-13 


(0x6) 





-12 


(0x6) 





-4 


N=7 


Ferro. 


(0.19,0.19,0.32x5) 


2 





(0.21,0.21,0.31x5) 


2 





(0.20,0.20,0.32x5) 


2 





AF1 i 


^-0.36,0.36,-0.33,-0.22,0.22,0.32,0) 





9 








(-0.32,0.32,-0.30,-0.20,0.20,0.30,0) 





8 


AF2 


(0,0,-0.36,-0.23,0.23,0.36,0) 





14 


(0,0,-0.29,-0.20,0.20,0.29,0) 





8 


(0,0,-0.32,-0.21,0.21,0.32,0) 





12 


Radial 


(0.27,0.27,0.18x5) 





22 


(0.24,0.24,0.12x5) 





14 


(0.24,0.24,0.17x5) 





20 


Para. 


(0x7) 





37 


(0x7) 





24 


(0x7) 





33 



The Pd6 cluster displays a behavior different from the 
rest, and therefore we discuss it separately. Futschek 
and coworkers 0] found that Pd_R was also ferromag- 
netic in contrast to Kumar et al.[18|, who predicted it 
to be paramagnetic. We have found that both states are 
nearly degenerate, with the paramagnetic solution being 
slightly more stable. Aditionally, we have been unable to 
find non-collinear or antiferromagnetic solutions for this 
cluster. 

In contrast, and independently of the pseudopotential 
or approximation used, the rest of the clusters show a 
rich variety of antiferromagnetic and non-collinear solu- 
tions. Most of these solutions, though not all, exist for 
all LDAl, LDA3 and GGA1. We have also found that, 
whenever they exist, the relative order of the different 
solutions is maintained, and the size of the atomic mo- 
ments is very similar. These facts strengthen our belief 
that Pd atomic clusters are much more insensitive to the 



pseudopotential and approximation employed than bulk 
Pd. It is also reassuring that most of the collinear so- 
lutions have been identified in previous calculations [l^j 
(e.g.: AF1 for Pd 4 and Pd 5 and AF2 for Pd 7 ). 

The non-collinear solutions found can be classified 
into those that release antiferromagnetic frustration and 
therefore have lower excitation energy than the AF solu- 
tion (NCI and NC2 in Pd 4 , shown in Table II and Fig. 
2), and radial or quasi-radial solutions, that resemble the 
hedgehogs found in low dimensional theories of classical 
or quantum antiferromagnetsl^. Hedgehogs in these 
theories do not release frustration but rather are excita- 
tions over the antiferromagnetic ground state. We also 
find that these radial states have a higher energy that the 
antiferromagnetic solution, and therefore do not release 
frustration. 

Notice that the antiferromagnetic and non-collinear so- 
lutions can be reached at temperatures of the order of 



4 



exchange and correlation potential. 



G-O AF1 
Q-H AF2 
O-O Radial 
A-A Para. 




0.95 1 1.05 
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FIG. 3: Excitation energy per atom of the magnetic solutions 
of the Pds cluster as a function of the average interatomic 
distance, using GGA1. 

room temperatures (25 meV). Therefore, any measure- 
ment of the magnetization performed at room temper- 
ature should find a thermal average of all those states, 
many of which have a tiny magnetic moment. It should 
not be surprising that such a measurement give a small 
net moment. 

We finally discuss the relationship between magnetism 
and equilibrium interatomic distances. We have found 
that these are essentially the same regardless of the mag- 
netic state for the largest clusters (n = 5 - 7), the small- 
est ones showing slight variations of less than 0.04 A, but 
only within the LDA solutions. We have additionally an- 
alyzed the relative stability of the different solutions as 
a function of the interatomic distance. To this aim, we 
plot the energy per atom of the low-lying excited states of 
the Pds cluster, relative to the ground state energy, as a 
function of an uniform volume expansion, obtained using 
GGA1. The figure shows that no crossover takes place, 
apart from the nearly-degenerate AF1 and AF2 solutions, 
that cross at an expansion of about 4%. Moreover, the 
relative energy differences arc essentially preserved and 
the local magnetic moments kept constant, except for the 
AF2 and radial solutions, where they slightly change (by 
about 10%). 

To summarize, we have studied the geometry and 
magnetic properties of the ground state and lowest ly- 
ing isomers of small palladium clusters Pdjv, with iV 
ranging from 3 to seven. Our results confirm that the 
ground state is indeed collinear or paramagnetic. We 
have found a rich variety of non-collinear low-lying iso- 
mers, some of which efficiently release frustration, while 
other (hedgehog- like solutions) do not. All these solu- 
tions should contribute to the room temperature mag- 
netic behavior of the clusters, probably rendering small 
measured magnetic moments. We have finally found that 
all these states are rather insensitive to the choice of the 
pscudopotential and to the approximation used for the 
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